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We study the fluctuations in luminosity distances due to gravitational lensing by large scale (> 35 
Mpc) structures, specifically voids and sheets. We use a simplified "Swiss cheese" model consisting 
of a ACDM Friedman-Robertson- Walker background in which a number of randomly distributed 
non-overlapping spherical regions are replaced by mass compensating comoving voids, each with a 
uniform density interior and a thin shell of matter on the surface. We compute the distribution of 
magnitude shifts using a variant of the method of Holz & Wald (1998), which includes the effect of 
lensing shear. The standard deviation of this distribution is ~ 0.027 magnitudes and the mean is 
~ 0.003 magnitudes for voids of radius 35 Mpc, sources at redshift Zs = 1.0, with the voids chosen 
so that 90% of the mass is on the shell today. The standard deviation varies from 0.005 to 0.06 
magnitudes as we vary the void size, source redshift, and fraction of mass on the shells today. If 
the shell walls are given a finite thickness of ~ 1 Mpc, the standard deviation is reduced to ~ 0.013 
magnitudes. This standard deviation due to voids is a factor ~ 3 smaller than that due to galaxy 
scale structures. We summarize our results in terms of a fitting formula that is accurate to ~ 20%, 
and also build a simplified analytic model that reproduces our results to within ~ 30%. Our model 
also allows us to explore the domain of validity of weak lensing theory for voids. We find that for 35 
Mpc voids, corrections to the dispersion due to lens-lens coupling are of order ~ 4%, and corrections 
to due shear are ~ 3%. Finally, we estimate the bias due to source-lens clustering in our model to 
be negligible. 



I. INTRODUCTION 

A. Background and Motivation 

A number of surveys are being planned to determine 
luminosity distances to various different astronomical 
sources, and to use them to constrain properties of the 
dark energy or modifications to gravity that drive the 
cosmic acceleration. It has long been recognized that 
perturbations to luminosity distances from weak gravi- 
tational lensing will be a source of error for these stud- 
ies, both statistical and systematic [THS]. For supernovae 
the lensing noise becomes significant only at high red- 
shifts [5], but for gravitational wave sources the lens- 
ing noise dominates over the intrinsic luminosity scatter 
[7|[S]. Theoretical predictions for the magnification prob- 
ability distribution can be folded into the data analysis 
of surveys to improve the results [9J, and in particular it 
is possible to exploit the known non-Gaussian nature of 
this distribution 10 . In addition, it is possible to treat 
the "lensing noise" in luminosity distances as a signal in 
its own right, which provides useful information |11| . (A 
tentative detection of this signal in supernovae data has 
been claimed in Ref. [E]-) For these reasons, it is use- 
ful to have a detailed understanding of the magnification 



probability distribution. 

There are a number of methods that have been used to 
study the effects of weak lensing on luminosity distances: 

• Weak lensing theory can be used to predict the 
variance of the magnification distribution from the 
matter power spectrum |13j . However, the accu- 
racy of this approach is limited and in particular it 
does not allow one to probe the non-Gaussian tails 
of the distribution^ . 

• One can use numerical ray tracing using the results 
of cosmological simulations of large scale struc- 
ture, such as the Millennium simulation ^15j and 
the Coyote Universe project [16], see, eg. Ref. [17]. 
This approach is highly accurate and is based on a 
realistic density distribution. However it requires 
substantial computational power and is also lim- 
ited in some other respects. The largest simula- 
tions to date are are confined to a cube of comov- 
ing size z 0.16, so only a limited range of source 
redshifts can be considered. Although the calcula- 
tions evolve large scale structure nonlinearly, it is 
impractical to get a continuous description of the 
evolution, which is needed for computing the per- 
turbations to light ray paths; only snapshots of the 



^ We note however that there is a proposal for an approximate 

•Electronic address: flanaKan@a stro.cornell.edu| "universal probability distribution" for magnifications that takes 

tElectronic address: nk236@cornell.edu' as input only the variance of the distribution as predicted by 

tElectronic address: ira@astro.cornoll.cdu weak lensing theory, and which would allow prediction of the 

^Electronic address: rav@kicp.uchicago.edu non-Gaussian tails |14| . 
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density distribution are available. Finally, because 
the calculations required to evolve the matter dis- 
tribution are formidable, it can be difficult to com- 
prehensively survey the space of the underlying pa- 
rameters of the model, such as the primordial per- 
turbation spectrum. 

• A third approach is to use simplified analytical 
models of the distribution of matter that allow 
rapid computation of the full probability distribu- 
tion of magnifications, see, eg., Refs. [2l l4l [T8l fT9] . 



In this paper we follow the third approach. We develop 
an idealized "Swiss cheese" model pi EHVEH] of large scale 
structure to study the effect of density inhomogeneities 
on luminosity distances. Our model is complementary to 
many of the existing models in that we focus on lensing 
produced by structures at the largest scales, voids and 
sheets, rather than that produced by individual galaxies 
and halos, the focus of many existing models. 



B. Our void model 

In "Swiss cheese" models [1 [lOHlS], the Universe 
contains a network of spherical, non- overlapping, mass- 
compensated voids. The voids are chosen to be mass 
compensated so that the potential perturbation vanishes 
outside each void. We idealize these models even fur- 
ther by assuming that each void consists of a central, 
uniformly underdense region surrounded by a zero thick- 
ness shell. Mass flows outward from the evacuated in- 
terior and is then trapped on the wall. Although it 
would be more realistic to consider voids with smooth 
density profiles, this very simplified model should cap- 
ture the essence of the effect of large scale density inho- 
mogeneities on luminosity distances. Since voids in the 
observable Universe tend to be surrounded by shells that 
are relatively thin compared to the size of their evacuated 
interiors, the idealization of zero thickness may not be a 
severe simplification, particularly because we expect that 
the main effect of inhomogeneities on the luminosity dis- 
tance depends only on the integral of the density contrast 
along the line of sight from the source to the observer. 
A key feature of our idealized models is that they can be 
evolved in time continuously and very simply. 

Within the context of this highly idealized class of 
models, we study the distribution of magnitude shifts 
relative to what would be found in a smooth cold dark 
matter (CDM) model of the Universe with a cosmologi- 
cal constant. A, for different void sizes and present day 
interior underdensities, and for a range of different source 
redshifts. Moreover, although we shall use a Newtonian 
description that is valid as long as the void radii are small 
enough compared with the Hubble length Hq^, the cal- 
culations can be made fully relativistic if desired. (We 
discuss some corrections that are higher order in HqR, 
where R is the void radius.) 



This paper is a follow-up to our earlier work [50] 
(henceforth VFW08), in which we considered the effect 
of a randomized set of voids with a single and rather large 
comoving radius, 350 Mpc, using a particular model for 
a smooth underdense interior inside a mass compensated 
shell. That study found that for a source with redshift 
Zs = 1.8, the mean magnitude shift relative to smooth 
flat, CDM for an ensemble of realizations of large scale 
voids was unimportant (-0.003), but the distribution of 
magnitude shifts was fairly broad, with a standard devi- 
ation of about 0.1. Here, we consider a wider range of 
redshifts and void sizes, and compute magnitude shifts 
relative to a more realistic ACDM background with mat- 
ter density today = 0.3 and dark energy density 
today r^A = 0.7. 



C. Predictions for lensing noise 

Our results for the standard deviation <7m of the mag- 
nitude shifts are summarized by the approximate fltting 
formula 



cr,„ w (0.027 ± 0.0007) 



R 



35 Mpc 



0.9 



(1.1) 

Here R is the comoving radius of the voids, Zg is the 
source redshift, and /o is the fraction of the total void 
mass in its shell today. The exponents are a = 0.51 ± 
0.03, /3 = 1.07±0.04, 7 = 1.34±0.05. This flt is accurate 
to - 20% for 35 Mpc < R< 350 Mpc, 0.01 < /o < 0.9, 
and 0.5 < Zg < 2.1. The mean magnitude shift is again 
unimportant, roughly a factor of ten smaller than the 
standard deviation (1.1 1. 



Our result (1.1) is computed in the limit of zero shell 
thickness. This idealization is not very realistic, since 
as we discuss in Sec. |III| below there is a logarithmic di- 
vergence in the variance of the lensing convergence in the 
zero thickness limit. This divergence arises from rays that 
pass very near to the void walls. The variance in the mag- 
nitude shift, however, is finite because of the nonlinear 
dependence of magnitude shift on lensing convergence; 
the divergence is cut off at lensing convergences of order 
unity. (The divergence can also get regulated by finite 
sampHng effects; see Sec. III|. To address this issue we 



also consider a more realistic model with void walls of 
some finite thickness Ar. We estimate in Sec. IIIIDl that 
for /o = 0.9, i? = 35 Mpc, and Zg ~ 1.0, the standard 
deviation in magnitude shift is 



0.013W1 -t-0.231n 



IMpc 
Ar 



(1.2) 



a factor of ~ 2 smaller than the thin-shell limit ( 1.1 1 for 
Ar = 1 Mpc. 



The rms magnitude shift (1.2) due to voids is a fac- 
tor of ^ 3 smaller than that computed from individual 
galaxies and halos in accord with expectations from 
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weak lensing theory using the power spectrum of den- 
sity perturbations (see Ref. and Appendix [A]). Thus 
lensing due to voids is subdominant but not neghgible. 

We also use our model to estimate the sizes of vari- 
ous nonlinear effects that go beyond hnear, weak-lensing 
theory. We estimate that for i? = 35 Mpc voids, the 
dispersion am is altered by ~ 4% by lens-lens coupling, 
by ~ 3% by shear. There are also large nonlinearities 
30% — 40%) in our model that arise from the nonlin- 
earity of void evolution. These results are qualitatively 
in agreement with some previous studies of nonlinear de- 
viations from weak lensing theory [2^26] . 

We also study the source-lens clustering effect fTT , the 
fact sources are more likely to be located in high density 
regions, which enhances the probability of a lens being lo- 
cated near the source. We estimate that the correspond- 
ing bias in the distribution of magnifications is negligible 
in our model. 



A. Newtonian model of a single void 

As discussed in the introduction, we will consider void 
radii R ranging from 35 Mpc to 350 Mpc, which are small 
compared to the Hubble length. Therefore we can use 
Newtonian gravity to describe each void; the correspond- 
ing error is of order (HqR)^ <C 1 which we ignore. 

We choose the background cosmology in which we 
place our voids to be an FRW Universe with matter frac- 
tion Hm and cosmological constant fraction 1 — ^Im- We 
denote by acx{t) the corresponding scale factor, which 
is normalized so that a^^lt) = 1 today. It satisfies the 
Friedman equation 



M 



1 - n 



M 



(2.1) 



where Ho is the Hubble parameter, which has the solu- 
tion 



D. Organization of this paper 



SHot^l — flu J -1-1 
= smh 



' 3/2' 
3/2 



(2.2) 



This paper is organized as follows. Section [Il| reviews 
our Swiss cheese void model. We discuss how the voids 
evolve in an FRW background and describe the model pa- 
rameters. Next, we describe how our void locations are 
randomized, by choosing impact parameters randomly as 
light rays exit one void and enter the next. Finally, we de- 
scribe our method of computing the magnification. Sec- 
tion |11I| describes our simple analytical model which re- 
produces the results of the simulations to within ~ 30%. 
It also describes a modification of our void model in which 
the shell walls are given a finite thickness, and gives the 
corresponding analytical results. Section |IV] gives the re- 
sults of our Monte Carlo simulations for the the probabil- 
ity distributions of magnifications, and discusses the de- 
pendence of the variance on the parameters of the model. 
In Section |Vj we study the source-lens clustering effect 
and the associated bias. Section IVll summarizes our re- 
sults and their implications. In Appendix |A] we discuss 
the power spectrum of our void model and the corre- 
sponding weak lensing prediction. Appendix [B] reviews 
the derivation of the method we use to compute the mag- 
nification distribution. Finally, Appendix [C] is a compar- 
ison of our results with other recent studies of lensing 
due to voids [H [THl [MM 123] . Our results are broadly 
consistent with these previous studies but our model is 
simpler in several respects. 



II. SIMPLE MODEL OF LENSING DUE TO 
VOIDS 

In this section we describe our simplified Swiss cheese 
model of large scale voids, and explain how we compute 
the distribution of magnifications in the model. 



Here oa = (flM / (1 — ^m))^ is the scale factor at which 
the cosmological constant starts to dominate. 

Our void model consists of a spherical region of con- 
stant comoving radius i?, with a uniform density interior 
surrounded by a thin shell. We assume that the void is 
mass compensated, so the total mass enclosed is the same 
as what it would be in FRW, namely 



M = 



2G ■ 



(2.3) 



We denote by f{t) the fraction of this mass in the thin 
shell, so that the mass in the interior is [1 — f{t)]M. The 
fractional density perturbation in comoving coordinates 
(5m(x, i) = Sp{x,t)/p is therefore 

J^(x, t) = -/(i)e(i? - r) + ^fit)R6ir - R), (2.4) 

where Q{x) is the function defined by 8(a;) = 1 for x > 
and e(a;) = for x < 0. 

The corresponding potential perturbation 0, in a New- 
tonian gauge in which the metric has the form 



(1 + 20) dt^ + al^{t) {1 - 24^) 



(2.5) 



is given by solving the Poisson equation V^t; 
37?^riA/'5m(x, t)/(2acx)- This gives 



(x,<) = 



4aox(t) 

The corresponding radial acceleration is 



{R^ -r^)e{R-r). (2.6) 



2aox(i)' 



-re(i?-r) 
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For each void, the potential will take the form (2.6) in a 
spherical polar coordinate system centered on that void, 
and the total potential is given by summing over the 
voids. The potential vanishes in between the voids. 

We next discuss how to compute the time evolution of 
the fraction f{t) of the void mass in the thin shell. We 
will work to leading, Newtonian order in {HqR)^, and we 
will also neglect the surface pressure that would arise in 
a relativistic calculation. The uniform interior behaves 
like a positive energy FRW cosmology. It has negative 
curvature, fc < 0, and a scale factor ain(i) that obeys the 
equation 



1 - n 



k 



M 



(2.7) 



since the cosmological constant is the same inside and 
outside the void but the matter density is not. We de- 
fine the positive constant — ~^mHq /k, the inverse 
of which is proportional to the density contrast at early 



times. The solution to Eq. (2.7) is 



M 



dx 



(2.8) 



1 



This solution assumes that ai„ — a^x = at t = 0, so that 
the interior and the exterior regions started expanding at 
the same time. Otherwise the deviations from FRW are 



l arge at early times. Eliminating t between Eqs. ( 2.2 1 and 
(2.8) gives the relationship between a 



and flex, which is 



sinh 



' 3/2' 
3/2 



dx 



(2.9) 



1 + x2 



Note that the above equations imply that > flex, as 
fc < 0. The density of the interior is equal to the mean 
FRW density times (aex/oin)^ < 1, and so the fraction of 
mass in the shell is 



m 



1 



(2.10) 



We numerically solve Eq. (2.9 1 to obtain acx/ain as a 
function of Oex/aA, ^^_a/ and oq. In the remainder of the 
paper, we will parameterize our void models in terms of 
the value today /o ~ /(^o) of the mass fraction f{t) in 
the shell. We will usually pick /o — 0.9. The parameter 
flo can be computed from /o and ^m- 



B. Algorithm for randomization of void placement 

We now discuss how we choose the number and loca- 
tions of voids in our model. In some previous studies [34]- 
l36] . the centers of all the voids encountered by a given 
ray were chosen to be collinear, so that the ray passed 
through the centers of all the voids. In these studies the 



lensing demagnification was large enough to successfully 
mimic the effects of dark energy. However, as discussed 
in VFW08, the large demagnification was an artifact of 
the non-randomness of the void locations, which is not 
in accord with observations of the distribution of voids 
[37Vl40j . In this paper, we use a more realistic void dis- 
tribution, which we compute according to the following 
procedure: 

1. Fix the comoving void size R. 

2. Fix the redshift of the source Zg. 

3. Place voids all along the ray from the source to the 
observer, lined up so that they are just touching. 
The source and the observer are placed in FRW 
regions. The distance from the source to the shell 
of the adjacent void is chosen to be a fixed small 
parameter, and the distance between the observer 
and the shell of the adjacent void then depends on 
the number of voids that can fit between the source 
and observer. 

4. Randomize impact parameters by shifting each void 
in a random direction perpendicular to the direc- 
tion of the light ray, so that the square of the 
impact parameter is uniformly distributed between 
and 

Note that with this algorithm, each ray spends some 
time in FRW regions between each pair of voids. An 
alternative procedure would that used by Holz & Wald 
121, in which after exiting a void, a ray immediately en- 
ters another void without traversing an FRW region. In 
this model the effective packing fraction of voids would 
be a factor ^ 2 or so higher than in our model, and the 
rms magnifications and demagnification would be corre- 
spondingly enhanced. 



C. Method of computing magnification along a ray 

We now turn to a description of the method we use to 
compute the magnification for a ray propagating through 
a Universe filled with randomly placed voids, as described 
in the last subsection. Our method is essentially a modifi- 
cation of the method introduced by Holz & Wald [Ij , and 
goes beyond weak-lensing theory. In this section we de- 
scribe the computational procedure; a derivation is given 
in Appendix [B| 

Starting from the perturbed FRW metric ( |2.5[ ), we con- 
sider an observer aX t — to (today) and x = 0, or equiv- 
alently at rj — rjQ, where t] = J dt/ac^it) is conformal 
time. We consider a source at x = = Xs'n., where 
n is a unit vector. The geodesic joining the source and 
observer in the background FRW geometry is 



{x) = (770 - X, nx) , 



(2.11) 



for < a; < Xj, where x is the comoving distance (or 
affine parameter with respect to the flat metric ds^ = 



5 



aex {ty^ ds"^ = -df]"^ + dx2). Following Holz & Wald g], 
we solve the geodesic deviation equation relative to this 
unperturbed ray in order to find the net magnification 
and shear. We do not include deflection of the central 
ray since the resulting corrections are relatively small; 
see Appendix [b] and Ref. [2]. 

We denote by fc = d/dx = —djj + n^di the past directed 
tangent vector to the ray. We also introduce a pair of 
spatial basis vectors ca, A — 1,2, so that ca and n are 
orthonormal with respect to ds^. We define the projected 
Riemann tensor 



(2.12) 



for A, B — 1,2 where Ra^ps is the Riemann tensor of the 
perturbed FRW metric without the (t)'^ factor: 



ds^ 



(1 + 2</>) dr/^ + (1 - 20) dx^ 



Next we consider the differential equation along the 
ray 



dx^ 



{x) = -n%{x)A''B{^): 



(2.14) 



where 'R-'^q {x) means the projected Riemann tensor eval- 
uated at = x" (x), and capital Roman indices are 
raised and lowered with 5ab- We solve the differential 



equation (2.14) subject to the initial conditions at the 
observer 



A^B (0) 



0, 



dA' 



dx 



(0) 



5\ 



(2.15) 



Finally the magnification along the ray, relative to the 
background FRW metric, is'^ 



|det^(a::s)| ' 



(2.16) 



where the right hand side is evaluated at the location 
X — Xs oi the source. Note that this quantity, the ra- 
tio between the perturbed and unperturbed angular di- 
ameter distances, is a conformal invariant, as we show 
explicitly in Appendix [B} 

The matrix A{xs)/Xs can be expressed as a product 
of an orthogonal matrix and a symmetric matrix with 
two real eigenvalues 1 — k± 7, where k is called the lens- 
ing convergence and 7 the shear. The magnification is 
therefore 



(1 



■ 7 



(2.17) 



This computational procedure is essentially the same 
as that used by Holz & Wald |2j, except that Holz & 
Wald work in the physical spacetime rather than the con- 
formally transformed spacetime, and at the end of the 
computation they compute the ratio between the quan- 
tity x'j. / idei^) evaluated in the perturbed spacetime and 
in the background spacetime. In our approach we do 
not need to compute a ratio, and furthermore the source 



(2.13) 

we have 



term in the differential equation (2.14) vanishes in FRW 
regions between the voids, which simplifies the compu- 
tation. See Appendix [B] for more details on the relation 
between the two approaches. 

We now turn to a discussion of the method we use to 
compute approximate solutions to the differential equa- 
tion (2.14). Consider a small segment of ray, from x — xi 
to X — X2 say. Since the differential equation is linear. 







.A^b{^2)_ 





j\,{x2,Xi) K^^{x2,Xi) 
L^C Xi) M^p (X2, Xi) 

A^eixi) 



A'^six,) 



(2.18) 



for some 2x2 matrices J, K, L, M which together form 
a 4 X 4 matrix. To linear order in TZab we have'^ 



S%- dx{x2-x)n^ci^)^ (2.19a) 

{X2 - Xi) 



xi 

A 
C 

X2 rx 



1^ n ~ 



dx I dx(x-xi)7e^c(x),(2.19b) 

rX2 

/ dxTZ\, (x) ) (2.19c) 

J Xi 



S%- / dx{x-xi)n^cix)- (2.19d) 



We evaluate these matrices for a transition through a sin- 



gle void, using the potential (2.6|, the metric (2.13) and 



the definition (2.12) of TZab- We neglect the time evo- 
lution of the potential during passage through the void; 
the corresponding corrections are suppressed by {HqR)^. 



^ In our Monte Carlo simulations we discard all cases where the 
determinant is negative, and so the absolute value sign in Eq. 

can be dropped. As explained in Ref. 0, this prescription 
yields the distribution of magnifications of primary images; it 
is not possible using the geodesic deviation equation method to 
compute the distribution of total luminosity of all the images of 



Holz & Wald [2] dro p all o f the integrals over the projected Rie- 



a source. 



This is valid to leading order in {HoR)^. We keep the extra terms 
m Eqs. ( [2l9| even though our formalism neglects other effects 
that also give fractional corrections of order (HqR)^ . The extra 
terms change am by a few percent. 
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This gives 



1 4 
4 1 



ix2 ~x,)s\ + -c'r{z) 



1 2 

2 1 



L\ = 2c-p(z) 



«1 

3c2 



+ 20275 (z) 



3c2 



3c2 



3c2 



(2.20a) 
(2.20b) 

(2.20c) 
, (2.20d) 



Here b is the impact parameter, c — \/ — 6^ , 
7'(z) 



-HqHm — 



flex (2) ' 



(2.21) 



and /(z) is defined by Eq. (2.10). In these equations x 
and z are evaluated at the center of the void. 

Our computational procedure can now be summarized 
as follows: 

1. Pick some source redshift Zg, void radius R, and 
fraction of void mass on the shell today /q. 

2. Choose void locations according to the prescription 
described in Sec lII Bl 

3. For each void, compute the 4x4 matrix that is 
formed by the matrices J, K, L and M from Eqs. 



(2.20). 



4. Perform a similarity transformation J — >■ U"""^ • J-U 
on each of the matrices J,K,L,M for some ran- 
domly chosen SO{2) matrix U, to randomize the 
direction of the vectorial impact parameter. 

5. Multiply together all the 4x4 matrices, and mul- 
tiply by the initial conditions (2.15), to evaluate 



6. Compute the magnification /i relative to FRW from 
Eq. (2.16), and then distance modulus shift Am 
from 



Am 



logio (m) 



2 In 10 



In 1(1 



7^1. 



(2.22a) 
(2.22b) 



randomly on a sphere at a distance corresponding 
to redshift Zg. The corrected distribution is [2j 

7'(Am;z,) = J\fp{Am;Zs)/fi 

= A/'p(Am;z,)102^™/^ (2.23) 

where jV is a normalization constant. 

D. Relation to weak lensing theory 

In weak lensing theory the matrix Ji(xs)/xs that de- 
scribes the defiections of the rays is presumed to be al- 
ways very close to the unit matrix, so the total integrated 
effect of the inhomogeneities on a given ray can be treated 
linearly. The solution to Eq. (2.14) in this approximation 



is given by Eq. (2.19b) with xi = 0, X2 = Xg, 



dx 



X{XS — ^) rJ^A 



(2.24) 



Taking the determinant and linearizing again, the con- 
tribution from shear vanishes and the magnification is 
fi = 1 A- 2k where the lensing convergence k is given by 
the standard formula 



dx^^^^<5,n(x). (2.25) 

X s^cy. \Z) 

Here 5-^-i{x) is the fractional over density, x is the comov- 
ing distance, Xs is comoving distance to the source, and 
flex {z) is the scale factor. Evaluating this for our void 
model gives 



where the sum is over the voids and 



Xi {^Xs 



*^s^cx(^z) 



-f{Zi)Ci 



3c2 



(2.26) 



(2.27) 



is the lensing convergence from the ith void. Here Zi and 
Xi are the redshift and comoving distance to the center 
of the ith void, Ci = y/ R"^ — tif and bi is the ith impact 
parameter. Our model goes beyond the weak lensing 
result (2.26) as it includes lens-lens couplings and shear. 



III. APPROXIMATE ANALYTICAL 
COMPUTATION OF MAGNIFICATION 
DISPERSION 

A. Overview 



7. Repeat steps 2 to 6 a large number of times to gen- 
erate the distribution p{Am; Zg) of distance modu- 
lus shifts Am for sources at redshift Zs, for a ran- 
domly chosen direction from the observer. 

8. Finally, we correct this distribution to obtain the 
observationally relevant quantity, the probability 
distribution of magnitude shifts for a source chosen 



In the previous section, we described a Monte Carlo 
procedure for computing the probability distribution 
7^(^171; Zs) of magnitude shifts Am for sources at red- 
shift Zs, for our Swiss cheese model of voids. We will be 
particularly interested in the mean 

(Am) = dAmAmV{Am;Zs) (3.1) 



7 



and variance 

cr2„ = y dAm (Am- (Am))^7'(Am;z,) (3.2) 

of this distribution. In subsequent sections of the paper 
we wiU describe the results of our Monte Carlo simula- 
tions and their implications. In this section, however, we 
will take a detour and describe a simple, approximate, 
analytic computation of the variance. The approxima- 
tion consists of using the weak lensing approximation to 
compute the total lensing convergence k (accurate to a 
few percent, see Sec. IV B), and then using an approx- 



nonlinear relation (2.22b) between k and the magnitude 
shift Am. We will see in Sec. |IV| below that this analytic 
approximation agrees with our Monte Carlo simulations 
to within ~ 30%. 



Neglecting shear, the relation (2.22b) reduces to 
5 



Am 



In 10 



In 11 



(3.3) 



where k is given by Eqs. (2.26) and (2.27). We will see 
shortly that the the variance of k diverges. This diver- 
gence is an artifact of our use of a distributional density 
profile for each void, with a (5-function on the void's sur- 
face, and can be removed by endowing each shell with 
some small finite thickness Ar (see Sec. Ill Dj below). The 
variance of Am, on the other hand, is finite, because of 
the nonlinear relation (3.3). We shall proceed by using 
the linearized version 



Am 



5 

"in 10 



(3.4) 



of Eq. (3.3), and by simply cutting off the divergent in- 
tegrals that arise, at k ~ 1, the regime where the nonlin- 



earity of the relation (3.3) becomes important 



B. Variance of magnitude shifts 



From Eq. (3.4) we find for the mean and variance of 



the magnitude shift 
(Am) : 



In 10 



(inio) + Oi-')] ■ (3.5) 



The averages are over the set of impact parameters {6,; : 
i € [1, j(a:s)]} in Eq. (2.27), where j{xs) is the number of 
voids out to the source at Xs- In computing the averages, 
it will prove convenient to define 



1 - b^R' 



(3.6) 



so that each qi is distributed uniformly between zero and 
one, since impact parameters arbitrarily close to the void 
boundary are permitted. In fact, a shortcoming of our 



model is the vanishing thickness of the void wall. We 
therefore introduce lower cutoffs Ci for each void*, that 
is, we restrict Qi to lie in the range Ci < qi < I. We will 
discuss below the origin and appropriate values of these 
cutoffs. 

With this assumption we obtain for the mean of the 
lensing convergence (2.27) of the ith void 





\ r- 1 " 




3V'7j 7= 








-Q) 



imate cutoff procedure to incorporate the effect of the where 



Xi{Xs - Xi)f{zi){l + Zi) 



(3.7) 



(3.8) 



The mean lensing convergence (3.7) is always negative, 
since C'i < 1; introducing the cutoff leads to a bias toward 
de-magnification. This is a shortcoming of the model, 
since for any mass-compensated perturbation (k^) = ^. 
Below, we shall ignore small corrections that are powers 
of Ci, and will take (kJ = for all i. 

By contrast the second moment {nf) diverges logarith- 
mically in the limit Ci — > 0: 



[H^VLMXaRWiY j dqi(^^ 
(HlVLMXsRwi) 



InC; 



1 

0{C,) 



(3.9a) 
(3.9b) 



This divergence is caused by rays that just graze the 5 
function shell of the void. 

Because k is a sum of k^, its mean is the sum of the 
individual means, but 



Ki) {Kj 



and therefore 

a^ = (.2)-(«)^=5^((.?)-(..)^). 



(3.10) 



(3.11) 



Combining this with Eqs. (3.4) and (3.9b) and dropping 



terms linear in Ci gives for the variance in magnitude 
shift 



2 \ ^ 2 



-InC,- 



(3.12) 



* These cutoffs will be used only for construction of our analytical 
model in this section; they are not used in Monte Carlo simula- 
tions in the remainder of the paper. 

^ To restore this feature we could either scale the contribution from 
the underdense core downward by a factor of Si = 1 -I- -I- Ci 
or scale the contribution from the overdense shell upward by the 
same factor Si. 
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where we have defined 



Co 



h^io ■ 



(3.13) 



We choose the cutofl^s Ci to correspond to ~ 1, as 

this gives 



discussed above; from Eqs. (2.271 and (3 



(3.14) 



The approximate analytic resuh given by Eqs. (3.8) and 
(|3.12|) - (|3.14|) is plotted in Fig. 



10 in Sec. IV B below 



It agrees with our Monte Carlo simulations to within 
~ 30%, which is reasonable given the crudeness of our 
analytic cutoff procedure. 



which is properly normalized for Ci > qi^c if we extend 
the range of Ci to infinity, thereby incurring an error 
- exp(-iV). 



We now average the expression (3.12 ) for the width cr„j, 



using the distribution (3.18) to average over the cutoffs 
Ci. The result is 



m,mcd 



If we define 

S{fo,Zs) = E' 



(3.19) 



3 f°° 

In TV - - - / dxe""" \n{x + iVgi,c) 

2 Jn 



C. Finite sampling effects 

In addition to computing the width cr^,j of the distribu- 
tion of magnitude shifts Am, we now compute a different 
quantity 

incd('^) which is, roughly speaking, the esti- 
mate of the width that one would obtain with N samples 
Amjj, I < a ^ N , drawn from the distribution. More 
precisely, this quantity is defined as follows. From the N 
samples we construct the estimator 



1 



N 



a=l 



1 



N 



N-l^ iV(iV-l) 

a=l ^ ' \a=i / 

(3.15) 

This quantity is itself a random variable with expected 
value ((T^„) = cr^- However for finite N the median value 
of the distribution of ct^j can be significantly different 
from cr^. We denote this median value by cr^mod(-^)- 
In the limit iV — oo we have am,med(,N) — am- We note 
that realistic supernovae surveys will have no more than 
lO'* supernovae. 

To estimate this median value, we use the fact that for 
each void i, finite sampling imposes a minimum value on 
qi of qi ~ 1/A'^ on average, which acts like a statistical 
cutoff in the integral (|3.9b|. This is in addition to the 



physical cutoff (3.14) discussed above, which we will de- 
note by qi^c from now on. For N samples qi^ai ^ ^ ct < N , 
the probability that all N samples are larger than a value 
Ci which is larger than qi^c is 



Po(< C,) 



1-C, 

1 - qi^c 



N 



(3.16) 



Differentiating once we find that the probability distri- 
bution of C, is 



dPoi< a) 



N{1 - C) 
(1 - 9j:,c) 



N-l 



N 



(3.17) 



For very large values of N and small q^.c an adequate 
approximation is 



P(a) «Afexp[-7V(a-g,,e)], 



(3.18) 



liNq^,, 



dxe~''\ii{x + Nq.^c), (3.20) 



then Eq. (3.19) becomes 



2 

'^m,mcd 



Sifo,z,)[\iiN 



wh{Nq.,c) 

i 

(3.21 



The result (3.21) was obtained by averaging over the 



cutoffs {Ci} using the probability distribution (3.18), and 
is an estimate of the median of the distribution of (t^. 
Of course the actual value of (T^„ computed from a Monte 
Carlo realization of N lines of sight, or obtained from 
N observations of magnifications, may differ from the 
result (3.21). We would li ke to also estimate the spread 
in values of ct^ . From Eq. (3.12), and taking the variance 



with respect to the distribution of cutoffs Ci, we find 



\ _ w'^Var(iVq,,e) 
[E.^f(E^ + 3/2)]^ 



(3.22) 



m.mcd 



where 



Var(7VQ,,,) = (lnC,)2- (Ina) 

dxe--[\n{x + Nq,^,)f - [l{Nq,^,)Y 

(3.23) 

Here the overbars denote an expectation value with re- 



spect to the probability distribution (3.18). The quantity 



(3.22) is a measure in the fractional spread in our esti- 



mate of the median, and should give a lower bound on 
the fractional spread in values of (t^ 



Two limits of Eqs. (3.21) and (3.22 1 are especially 
simple. First, for Nqi^c ^ 1, we have ^{Nqi^c) ~ 
7£; = 0.5772 . . ., the Euler-Mascheroni constant, and also 
Var(iVg.,^c) ~ 1-645 and -InC; sa In + 7_e. This gives 



'^m,med 
^'^m.mcd 



a;^5(/o,z,)( In iV--+7i;), (3.24a) 



1.645 



cr, 



m,mcd 



iVvoid(lniV + 7B-3/2)' 



(3.24b) 
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where A^void = Xs/{2R) is the number of voids and we 
have used the crude approximation Wi = constant in 
the second equation. Second, for Nqi^c ^ 1, we have 
l{Nq^,c) ~ -lniVg,,c, Yai{Nq,^^) ^' l/{Nq,^^)^ , and 
InCi ~ lnq,;,c, and so we obtain 



m,mcd 



-S{fo,Zs) ~^Wi\nq,^c 



Aa. 



2 

m,incd 



m,incd 



N 



(3.25a) 
(3.25b) 



The second case (3.25a) coincides with the N- 
independent width (3.12) - (3.14) computed earher. We 



see that the resuhs are dictated by a competition between 
statistical and physical cutoffs via the dimensionless pa- 
rameter Nq.i^c- 

As discussed above, our simulations are effectively cut 
off at ^ 1; this implies a physical cutoff 



Wi 

2.2 X 10 



0.23 



V35Mpc 



(4mOt3.26) 



Here we have scaled the factor HmHqXs to its value at 
r^Af = 0.3, Zs = 1.0, the quantity ho.7 is given by Hq ~ 
70ft-o.7 kms~^Mpc~^, and we note that Awi < f{zi){l + 
Zi). From the estimate (3.26) we expect the Nqi^c ^ 1 
limit to apply for N < 10°. In this case the cutoff is 
purely statistic al and the physical cutoff is unimportant. 
The prediction (3.24a) for cr„i,niod for N = and Zs = I 
is shown in Fig. [1| together with results from our Monte 
Carlo simulations, which are described in Sec. IV below. 
The plot shows good agreement between the model and 
the simulations. 

For this CcLS6, 3- lower bound on the fractional spread in 
the values of a'^ around its median value is given by Eq. 



(3.24b). That is, in any given simulation or observational 



survey with N light sources, the scatter of values about 
the expected will be at least this large. For example, with 



1 and i? = 35 Mpc, the implied spread 



iV = 10*, z 

is > 6%. In this regime where the cutoff is primarily 
statistical, the range of likely values of am is substantial, 
and only decreases logarithmically with increasing TV. 
When TV > 10^, we move into the Nqi^c ^ 1 regime 



where Eqs. (3.25) apply. The results in this regime were 
discussed in Sec. IIIB| above, and are plotted in Fig. [TO] 



in Sec. IV B below. Equation (3.25b) indicates that the 
spread scales as 1/A^ in this regime. However this es- 
timate is only a lower bound for the spread in va lues 
of (T^, as discussed above. In fact, from Eq. (3.15) the 
standard deviation of ct^ can be computed in terms of N 
and of the second and fourth moments of Am; it scales 
like 1/Vn as TV — ?> oo. In any case, the spread decreases 
more rapidly as TV increases after the transition to the 
large TV regime. We will see in Sec. |IV| below that this 
prediction agrees well with our Monte Carlo simulations. 



0.010 
0.005 



0.001 
5 X IQ-* 




FIG. 1: The green line is our analytic model ( |3.24a[ ) of the 
median width of the distribution of magnitude shifts Am, for 
TV — 10* samples, source redshift Zg = 1.0 and void radius 
R — 35 Mpc, as a function of the fraction of mass /o on the 
void shells today. The data points are from our Monte Carlo 
simulations with the same parameter values, described in Sec. 
IIV] below. 



D. Extension of void model to incorporate finite 
shell thickness 



In this subsection we consider a modification of our 
void model, in which the void wall is given a finite co- 
moving thickness Ar^ that acts as a physical cutoff in 
the divergent integral (3.9a). The corresponding value 



of the cutoff parameter qi ^ is qi.c = 2Ari/i?, from Eq. 
(3.6). The value of wall thickness that corresponds to 
the cutoff (|3.26[) is thus Ar, ~ 3 pc(i?/35 Mpc)^, which 



is much smaller than the expected void wall thicknesses 
~ Mpc in large scale structure. Thus, our thin-shell void 
model is somewhat unrealistic; the results are modified 
(albeit only logarithmically) once the wall thickness ex- 
ceeds ~ pc scales. This motivates modifying the model 
to incorporate a finite wall thickness. 

Consider next how the wall thickness evolves with red- 
shift. At very early times, when the perturbation is in 
the linear regime, it maintains its shape in comoving co- 
ordinates, so the cutoff scale is some fixed fraction of 
R. Once the perturbation becomes nonlinear, the shell 
thickness should freeze out in physical extent, implying 
a comoving size oc 1/a. Thus, a suitable model for the 
redshift dependence of the cutoff would be 



qc{a) = eoW{a/aa), 



(3.27) 



where W{x) is a function with W{x) — > 1 for a; ^ 1 and 
W{x) — > Kq/x for X ^ I. Here ao(/o) is the scale factor 
when the perturbation ceases to be linear, and K and eg 
are constants that may also depend on /q. Very roughly, 
we expect qda) ^ 0.1 so Nqi^c ^ 1 as long as TV > 10, 
so that Eq. ( 3.25| ) will apply. 

Suppose now that for a restricted range of source red- 
shifts it suffices to take the fractional shell wall thickness 
Eg = Aril R in comoving coordinates to be the same for 
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FIG. 2: The magnitude shift Am as a function of source 
redshift Zg for a single run, for voids of radius i? = 35 Mpc, 
fraction of mass on the shell today /o = 0.9, in a ACDM 
cosmology with Q,m = 0.3. 



all shells. Then from Eq. \i.l2\ we get^ 

al^ = alS{h, zs) [- In + ln(2) + 0(e, In e,)] . (3.29) 



Equation (3.29) has the same form as Eq. (3.24a), but 
since Nts ^ 1, the implied am is smaller. For example, 
evaluating this expression for /o = 0.9, Zg = 1.0 and 
i? = 35 Mpc gives 



0.013a 



0.23 In 



lMpc\ 



(3.30) 



where Ar — tsR is the wall thickness. 

The logarithmic divergence of ct^ will also be regulated 
by treating the shell as composed of fragments that rep- 
resent local density enhancements such as galaxy clusters 
and superclusters for purposes of computing the magni- 
fication of passing light beams. We shall examine this 
further refinement of our model elsewhere. 



IV. RESULTS OF MONTE CARLO 
SIMULATIONS FOR MAGNIFICATION 
DISTRIBUTIONS 

We now turn to describing the results of our Monte 
Carlo simulations based on the algorithm described in 



Eq. ( |3.29| differs from Eq. ( |3.12| in that the -3/2 has been re- 
placed by In 2. This shghtly more accurate version of the equa- 
tion is derived as follows. Instead of using the cutoff procedure 
embodied in Eq. ( |3.9a[ l, we use a regulated density profile of 
the form 5^(r) = -fe{Ri - r) + ae{R - r)e(r - Ri) where 
Ri = — ts) and a = /[(I — es)~^ — The variance in 

the lensing convergence can then be computed from 



xrr In 



r^(^s s^i)"^ 

's '^ex (^■^i)'^ R"^ 



rR rR 
/ 

JO Jo 



(3.28) 
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FIG. 3: The probability distribution of magnitude shifts Am 
for a simulation in a ACDM cosmology with Qm ~ 0.3, with 
sources at redshift = 1, comoving voids radius R — 35 
Mpc, and fraction of void mass on the shell today fo = 0.9. 



Sec. |TT] In the remainder of this paper, unless otherwise 
specified, we will adopt the fiducial parameter values of 
matter fraction f2j\/ — 0.3, source redshift Zg — 1.0, void 
size i? = 35 Mpc, and fraction of void mass on shell today 
fo — 0.9. Our choice of void size is motivated by the fact 
that observed void sizes [57Hi5] range from a typical size 
of ~ 10 Mpc to an upper limit of ^ 100 Mpc. For this 
fiducial case, we show in Fig.[2]the distance modulus shift 
Am as a function of redshift Zg for a single realization of 
the void distribution. The values jump discontinuously 
after each void, and illustrate the stochastic nature of the 
lensing process. 

Next, we repeat this computation some large number 
N of times in order to generate the distribution of mod- 
ulus shifts Am. In the rest of the paper we will focus in 
particular on the mean (Am) and standard deviation cr„ 
of this distribution, and also on the estimator drn{N) of 
the standard deviation that one obtains at finite N ^ given 
by Eq. ( 3.15^ , which satisfies am{N) — )• am as N ^ oo. 



The distribution for the fiducial case for = 2 x 10 is 
shown in Fig. [3] For this case the standard deviation is 
am = 0.03135 ± 0.0003 and the mean is (Am) = 0.004 ± 
0.001 (where the error is estimated based on dividing the 
data into 200 groups of 10000 runs). Our result for the 
standard deviation agrees to within ~ 30% with that of a 
different Swiss cheese void model by Brouzakis, Tetradis 
and Tzavara ^29j : see Fig. 5 of that paper which applies 
to i? = 40 Mpc voids at = 1. It also agrees to within 
a factor ^ 2 with the predictions of weak lensing theory 
using an approximate power spectrum for our void model, 
as discussed in Appendix \K\ 

Figure |4] shows how our estimated standard deviation 
am{N) varies with number of runs N . The quantity plot- 
ted is logiQ |(7„j/(T„j — 1|, where cr„i = 0.03135 is an esti- 
mate of the A^ — > oo limit, here taken from our largest 
run with A^ = 10^. This plot exhibits several interesting 
features that are in good agreement with the analytical 



from Eq. i 2.25 I 



model described in Sec. Ill First, in the low A^ regime at 
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FIG. 4: The estimator a™, of the standard deviation of the 
distribution of magnitude shifts Am, as a function of number 
A'' of runs, for sources at redshift Zs = 1, comoving voids 
radius i? = 35 Mpc, and fraction of void mass on the sheU 
today fo = 0.9. The plotted quantity is logj^Q |(Tm/(Jm, — 1|, 



where Um = 0.03135 is an estimate of the A'' 
taken from our largest run with A'^ = 10*'. 



oo limit, here 
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say N ~ 10"*, the values of am(iV) differ systematically 
from the asymptotic value by a few tens of percent, re- 
flecting the difference between CTm.med and <Tm- Second, 
there is a somewhat smaller scatter in this regime, of 
~ 5%, in agreement with the prediction (3.24b I. Third, 
there is a transition to a different behavior at ^ 3x lO'^, 
after which both the scatter and systematic deviation 
from the asymptotic value are much smaller. 

In the rest of this paper, we will use the value N = 10^ 
unless otherwise specified. From Fig. [4] this corresponds 
to an accuracy of ~ 1 percent. 

We show in Fig. [sjthe mean (Am) of the distribution 
as a function of source redshift Zg, for i? = 35 Mpc and 
N = 2x 10^. The errors shown are estimated by dividing 
the data into 200 groups of 10000 runs. The effect of the 
nonzero mean on cosmological studies cannot be reduced 
by using a large number of supernovae, unlike the effect 
of the dispersion am- However, the mean (Am) ~ 0.003 
magnitudes shown in Fig. [5] is too small to impact cos- 
mological studies in the foreseeable future. 

In Figs. |6] [Tj and [8] we show the probability distri- 
butions of magnitude shifts Am for some other cases: 
source redshifts of Zg — 1.1,1.6 and 2.1, and void radii 
of i? = 35 Mpc, 100 Mpc, and 350 Mpc. We now turn 
to a discussion of the dependence of our results on these 
parameters, as well as on the fraction of mass in the shell 
today fo. 



A. Dependence on void size 

In Fig. [9] we show the standard deviation am of the 
magnitude shift as a function of void size i?, for three 
different redshifts, Zg = 1.1, 1.6, 2.1. To a good approxi- 
mation the standard deviation grows as the square root 
of the void size, am oc \/R. We can understand this 



FIG. 5: [Top] The mean (Am) of the distribution of magni- 
tude shifts Am as a function of source redshift Zs, for voids 
of radius i? = 35 Mpc with fraction of mass on the shell to- 
day fo — 0.9, for A^ = 10^ samples. [Bottom] The same for 
R = 100 Mpc. 



scaling by making some order of magnitude estimates. 

In making these estimates, we consider two different 
classes of rays. Consider first rays that never come very 
close to the shell of any of the voids, i.e. we exclude the 
case b — R R, where b is the impact parameter. The 
potential perturbation Ac/) for passage through a void is 
of order A0 ~ fR^H^, where / is the fraction of void 
mass in the shell (or equivalently the fractional density 
perturbation in the void interior). The contribution to 
the lensing convergence from this void is then of order 
K ~ A(l)/ (HqR) ~ fHoR. Next, the trajectory of rays is 
a random walk, so the net lensing convergence is the rms 
convergence for a single void multiplied by the square 
root of the number ^ 1/{HqR) of voids. Thus the con- 
tribution to the rms magnitude shift from this class of 
rays is of order 



fVH^- 



(4.1) 



Consider next rays which just graze the shell of at least 
one of the voids. These grazing rays are subject to large 
deflections, because of the (5-function in density on the 
surface of the void. The large deflections cause cause the 
second moment (k^) of the lensing convergence to di- 
verge, as discussed in Sec. IIIB However, the standard 
deviation of the magnitude shift Am is still finite, be- 
cause of the logarithmic relation (3.3) between Am and 

K. 

For estimating the effect of these grazing rays, we ne- 
glect shear. The convergence k of the grazed void will 
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FIG. 6: The probability distributions of magnitude shifts Am 
for simulations with sources at redshifts of Zs — 1.1 (top), 
Zs — 1.6 (middle) and Za = 2.1 (bottom), for comoving voids 
of radius R =35 Mpc with 90% of the void mass on the shell 
today. 



FIG. 7: The probability distributions of magnitude shifts Am, 
at source redshifts Zs of 1.1 (top), 1.6 (middle) and 2.1 (bot- 
tom), as in Fig.[6]except with comoving void radius oi R = 100 
Mpc. 



B. Dependence on fraction of void mass on the 
shell 



be of order unity or larger if the impact pa ramet er b is 
b = R{1- e), where e - pR^H^, from Eq. \2.27\ . This 
will occur with probability ~ e. The contribution of these 
rays to ((Am)^) oc ([ln(l — will be of order e times 
the number ^ 1/{HqR) of voids, or a,n ^ f\/HoR, the 



same as the result (4.1 1 for the non-grazing rays. 



These considerations show that both the underdense 
void and the mass-compensating shell make substantial, 
comparably large contributions to Om ■ This suggests that 
it may be important to refine the shell model to include 
its fragmentation into localized overdensities representing 
galaxy clusters and galaxies, as discussed in Sec. |IIID| 
above. 



In this subsection we discuss the dependence of the 
magnification distribution on the fraction /o of void mass 
on the shell today, or, equivalently, on the fractional over- 



density 5p/p, cf. Eq. (2.4) above. Figure 10 shows the 



results of our simulations for (7 777, RS 3t function of /o for 



N = W", together with a fit of the form (4.2) 



.(/o)=«/o-)-/?/2 



(4.2) 



for some constants a and /3. We find that a — 0.025 ± 
0.006 and /3 = 0.0085 ± 0.0064. Thus, the data show a 
statistically significant deviation from linear behavior, of 
the order of - 30 - 40%. 

We now discuss the various sources of nonlinearity that 
arise in the computation. We will consider three different 
types of effects. 

First, in weak lensing theory, the magnification is a 
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FIG. 8: The probability distributions of magnitude shifts Am, 
at source redshifts Zs of 1.1 (top), 1.6 (middle) and 2.1 (bot- 
tom), as in Fig.[6]except with comoving void radius of iZ = 350 
Mpc. 



linear function of the density perturbation. Our com- 
putation includes some nonlinear effects that go beyond 
weak lensing theory, specifically lens-lens coupling (the 
fact that the deflection due to one lens modifies the de- 
flection caused by subsequent lenses) and shear (the effect 
of the non-trace components of the matrices T?,"^^ and 
g). To explore the magnitude of these effects, we per- 
formed Monte Carlo simulations where we compute the 
lensing convergence for each vo id and add these to obtain 
the total lensing convergence (2.261, and then compute 



Am from k using the exact nonlinear relation (2.22b) 
for zero shear. The resulting value of a 



1, i? = 35 Mpc, N = W is cr„ 



for /o = 0.9, 
0.0292, about 7% 



smaller than the value am — 0-0314 obtained by multi- 
plying the 4x4 matrices. Thus, there is a ~ 7% change 
from lens-lens coupling and shear. For R — 100 Mpc, the 
change due to lens-lens coupling and shear is ^ 10%. We 
also performed simulations where we kept just the trace 
part of the matrix TZab, in order to exclude the effects 




FIG. 9: The standard deviation am of the distribution of 
distance modulus shifts Am as a function of void radius 7?, 
computed using A'^ = 10^ runs for each point. The bottom 
line (stars) is for sources at = 1.1, the middle line (squares) 
is Zs — 1.6, and the top line (diamonds) is Zs = 2.1. Void radii 
range from 35 to 350 Mpc and the fraction of void mass on 
the shell today is /o = 0.9. The lines are fits of the form 




FIG. 10: The standard deviation am as a function of the 
fraction fo of the void mass on the shell today, for void radii 
of i? = 35 Mpc and source redshift of = 1, computed using 
N — runs for each point. The dashed blue curve is a 
fit of the form am = q/o + Pfo- This plot shows that there 
are nonlinearities present at the level of ~ 30 — 40%. The 
solid green curve is the analytic model (3.121 - (3.141, which 
is accurate to ~ 30%. 



of shear, but included lens-lens couplings by computing 
4x4 matrices for each void and multiplying all these 
matrices. In this case the deviations of (t,„ from the full 
simulations are ^ 3% for fo — 0.9, Zg = I, R — 35 Mpc 
and ~ 6% for R = 100 Mpc. Thus, corrections due to 
shear are of this order. 

These nonlinearities due to lens-lens coupling and 
shear are significantly smaller than the nonlinearity 
shown in Fig. |10[ Thus other sources of nonlinearity 
must dominate. For the remainder of this subsection we 
will neglect lens-lens coupling and shear, to simplify the 
discussion. 

A second type of nonlinearity present in our computa- 
tions is the fact that the void mass fraction f{z) at some 
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FIG. 11: The factor h{z, fo) by which nonlinear evolution cor- 
rects the growth function D+ (z) of linear perturbation theory, 
for our void model. The upper curve is for /o = 0.9 and the 
lower curve is for fo — 0.5. 



redshift z depends nonlinearly on its value /q = /(O) 
today, due to nonlinearity in the void evolution. There- 
fore, even if we make the weak-lensing approximation of 
a linear dependence of the magnification on the density 
perturbation f{z), the magnification will still be a non- 
linear function of /q. We can parameterize this nonlinear 
evolution effect by writing 



fiz;fo) = .foD+{z)h{zJo), 



(4.3) 



where (z) is the growth function of linear perturbation 
theory, normalized so that Z?+(0) = 1, and the function 
h{z, /o) incorporates the nonlinearity. This function sat- 
isfies h{z, fo) — >■ 1 as /o — and also as z — 0, and can 
be computed using the results of Sec. 11 A| above. Figure 
[TT] plots this function for fo = 0.5 and fo — 0.9, and 
shows that the nonlinearities in the evolution are signifi- 
cant. 

This nonlinear evolution effect is the dominant source 
of nonlinearity in our simulations. To illustrate this, we 
define, for a given source redshift z^, the parameter 



/„id = /(^./2,,/o). 



(4.4) 



In other words, /mid is the fraction of void mass on the 
shell for voids halfway to the source, the distance where 
most of the lensing occurs. We can use /mid instead 
of fo as a parameter to describe our voids. With this 
choice of parameterization, the nonlinear evolution ef- 
fect is significantly reduced. This is illustrated in Fig. 
[12] which shows the same data as in Fig. [lO) but as a 
function of /mid rather than fo- The best fit parame- 
ters in the quadratic fit = a/mid + /^/mid ^'^^ 
a = 0.032 ± 0.005, /3 = 0.0016 ± 0.0057, showing that 
there is no statistically significant nonlinearity. 

A third type of nonlinearity in our simulations arises 
from the nonlinear relation between the lensing conver- 
gence K and the magnitude shift Am. This effect should 
be present in our data but is quite small. If we neglect 
lens-lens coupling, shear, and the nonlinear evolution ef- 
fect, then we expect logarithmic terms in the relation 
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FIG. 12: The standard deviation am as a function of the 
fraction /mid of the void mass on the shell for voids halfway to 
the source, for void radii of 7? = 35 Mpc and source redshift of 
Zs = 1. The solid line is a fit of the form am = a/mid + Pf^id- 
For this choice of parameterization there is no statistically 
significant nonlinearity detectable in the data. 



between am and /o, of the form 

al^af^+Pfl\nfo + ..., (4.5) 

where a and /3 are constants which are independent of 
fo- This follows from the analysis of Sec. IIIB above, 



where the logarithmic divergence in the variance is cutoff 
at K ~ 1; see Eqs. (3.12| and (3.14). However our data 



show that the logarithmic terms in Eq. (4.5) are quite 
small. 

Next, we discuss the effects of allowing a distribution 
of values of void mass fraction on the shell fo in our 
simulations, rather than having a fixed value. We per- 
formed simulations where we pick a value of / for for each 
void crossing according to the following prescription. We 
choose a random values for I/oq from a Gaussian distri- 
bution with a mean of 8 and a variance of 30, truncated 
to lie in the range that corresponds to < / < 1. Figure 
|13| compares the probability distributions for magnitude 
shifts with and without variations in /. Treating / as a 
random variable increases the standard deviation Um by 
-3%. 



C. Dependence on source redshift 



Figure 14 shows the standard deviation cr„i of the mag- 
nitude shift distribution as a function of source redshift 
Zg, for three different void sizes. The standard deviation 
increases with redshift faster than Zg. This increase is 
due in part to the increasing number of voids but there 
are additional factors. 

To understand the redshift dependence analytically we 
use the expression for the dispersion in lensing conver- 
gence from weak lensing theory, given by Eq. \A2\ in 
Appendix [A] The matter power spectrum A(fc, z)^ for 
our void model is proportional to fiz)"^, so we obtain 
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FIG. 13: A comparison of the probability distributions of 
magnitude shifts Am in two different cases: fraction of mass 
on the shell today fixed at fo — 0.9 (circles), and fo drawn 
from a distribution as described in the text (stars). In both 
cases void radius is i? = 35 and source redshift is Zs — 1.0. 
The spread in the shell surface densities gives rise to a wider 
distribution of magnitude shifts, by about ~ 3%. 
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FIG. 14: The standard deviation am as a function of source 
redshift Zs, computed using TV — lO" runs, for voids of radii 
R = 35 Mpc (red, crossed circles), 70 Mpc (green, squares), 
and 105 Mpc (blue, circles). The lines are fits proportional to 
the analytic estimate (4.6 1. 



that 



(k^) oc / dxw{x, Xs)'^ f{z)'^ 



(4.6) 



where w{x, X s) = ( 1 + z)Hqx{xs — x)/xs and f{z) is de- 
fined by Eq. ( |2.10[ ). In the range of redshifts 0.5 < Zg < 
1.5 this redshift dependence is approximately a power 
law, proportional to z^ '^^, to within ~ 5% percent''. This 
redshift dependence agrees with the results of our simu- 
lations shown in Fig. 14 to within ^ 10%. 



^ The asymptotic behavior at large Zs is that the expression ||4 
increases linearly in Zs- 



We complete this part of the analysis by giving a three 
parameter fit for the standard deviation function 
of void radius R, fraction of void mass on the shell today 
/o, and source redshift Zg. The result is 



am ~ (0.027 ±0.0007) 



R 



A 

0.9 



iro)' 



35 Mpc, 

(4.7) 

where the parameters are a = 0.51 ± 0.03, f3 = 1.07 ± 
0.04, 7 = 1.34 ± 0.05. This fit is accurate to - 20% for 
35 Mpc <R< 350 Mpc, 0.01 < fo< 0.9, 0.5 < z^ < 2.1. 



V. BIAS DUE TO SOURCES OCCURRING 
PREFERENTIALLY IN HIGH DENSITY 
REGIONS 

For sources which are randomly distributed in space, 
it is known that the total expected apparent luminosity 
of a source, including all primary and secondary images, 
must agree with that of the background FRW model [2]. 
Hence, in situations where the probability of caustics can 
be neglected, the probability distribution (2.23) of mag- 



nifications /X must be unbiased. Biases arise in our com- 
putations because of caustic effects, and also because we 
study the probability distribution of the magnitude shift 
Am, which is a nonlinear function of /i, cf. Eq. (2.22a). 



However, there is an additional fundamental source of 
bias which arises from the fact that sources are not ran- 
domly distributed in space, and instead preferentially oc- 
cur in high density regions, where they are more likely to 
be close to a lens. This is the source-lens clustering effect 
|27) . In this section, we make an analytical estimate of 
the bias 5m of the distribution of magnitude shifts that 
is due to source-lens clustering in our void model. 

In our computations so far in this paper, we have 
placed the source outside the voids, in the FRW regions. 
However, in reality most matter is concentrated on the 
edges of voids, and so sources are more likely to be on 
the void edges. If we demand that sources always be lo- 
cated on void edges, then the mean of the distribution is 
shifted by an amount (see derivation below) 



6m = 



1 



3 In (10) 



{l + z,)H^R^nMfs 



(5.1) 



Here Zg is the redshift of the source and fg — f{zs) is 
the fraction of mass on the shell for voids at the source 
redshift. Evaluating this estimate for = 0.3, Zg = 
1.0, i? = 35 Mpc, /o = 0.9 gives (5m ~ 5 x 10-^ and 
6m ^ 5 X lO^"' for R = 350 Mpc. These biases are below 
the accuracy of upcoming cosmology surveys. 

Turn now to the derivation of the formula (5.1). We 



start from the standard formula (2.25) for the lensing 



convergence in weak lensing theory. We consider just the 
contribution to k from the last void. In the integral, over 
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this void, we approximate the factors x and l/oox {z) as 
constants. Writing = Xg — x we obtain 



«^lastvoid — {1 + Zs)^m [ 

^ J last 1 



i]S„,{-K,t)dri. (5.2) 



We also neglect the time dependence of S^{x, t) for inte- 
grating over the last void. 

We now consider two different models for randomizing 
the relative displacement between the center of the last 
void and the source. We denote by b the transverse dis- 
placement of the void center from the line of sight, as 
before, and denote by rjy the distance from the void cen- 
ter to the plane through the source perpendicular to the 
line of sight. 

In our first model, we assume b and rjy are randomly 
distributed, proportional to bdrj^db, with < rjy < R 
and < b < R. Computing the integral (5.2) for our 



void model (2.4) gives 



^last void 



Hq^Im (1 



a 



(5.3) 

{r)v +a) Vv <a 



3q 



where a = \/R'^ — b'^. Now averaging over b 
and rjy gives the expected value of (^last void) = 
{l + Zs)H§R''nMfs/15. 

In the second model, we assume that b and 7].^ are 
correlated so that the source is always on the surface of 
the void. The average of Kiast void (b, 7]^) in this model is 



x'^last void/ 



i sinO K (RsinO, Rcosd) dO, 
Jo 



(5.4) 



which using the formula (5.3) gives zero. Subtracting the 



means of the two models gives an estimate of the bias, 
and multiplying the result by 5/ In 10 to convert from 5k 
to Sm gives the formula (5.1 ). 



due to 35 Mpc voids for sources at Zg = I is am = 0.031, 
which is ~ 2 — 3 times smaller than that due to galaxy 
clusters (see Appendix [A| below) . Also the mean magni- 
tude shift due to voids is of order 5m ^ 0.003±0.001. We 
also studied the bias that arises from the source-lens clus- 
tering effect, and estimated that the contribution from 
voids to this bias is quite small, of order 5m ~ 5 x 10~^. 
Refining our model by giving each void shell a finite thick- 
ness of ~ 1 Mpc reduces the dispersion dm by a factor 
- 2. 

We used our model to estimate the sizes of various non- 
linear effects that go beyond linear, weak-lensing theory. 
We estimate that for i? = 35 Mpc the dispersion am is al- 
tered by ~ 4% by lens-lens coupling, by ~ 3% by shear. 
For 100 Mpc voids these numbers become 3% and 6% 
respectively. 

Our simple and easily tunable model for void lensing 
can be used as a starting point to study more complicated 
effects. For example, one can use various algorithms to 
generate realizations of distributions of non-overlapping 
spheres in three dimensional space. Given such a realiza- 
tion one could use the algorithm of this paper to study 
correlations between magnifications along rays with small 
angular separations, which would be relevant to future 
pencil beam surveys [46]. Finally, our model is comple- 
mentary to other simplified lensing models in the litera- 
ture that focus on lensing due to halos but neglect larger 
scale structures, for example the model of Refs. [HI US]. 
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VI. CONCLUSIONS 

In this paper, we presented a simple model to study the 
effects of voids on distance modulus shifts due to gravi- 
tational lensing. A number of future surveys will gather 
data on luminosity distances to various different astro- 
nomical sources, to use them to constrain properties of 
the source of cosmic acceleration. The accuracy of the re- 
sulting constraints will be degraded somewhat by lensing 
due to nonlinear large scale structures. We studied this 
effect by considering a ACDM Swiss cheese cosmology 
with mass compensating, randomly located voids with 
uniform interiors surrounded by thin shells. 

We used an algorithm to compute the probability dis- 
tributions of distance modulus shifts similar to that of 
Holz & Wald |2J. The rms magnitude shift due to gravi- 
tational lensing of voids is fairly small; the dispersion am 



Appendix A: Comparison with weak lensing theory 

In this appendix we show that our results agree moder- 
ately well with the predictions of weak lensing theory, by 
computing an approximate matter power spectrum for 
our void model. We also obtain an independent estimate 
of the lensing due to voids by using the power spectrum 
of the Millennium simulation J^J. 

It is somewhat complicated to compute an exact power 
spectrum for our distribution of voids. As a simple 
model, we choose a two-void probability distribution 
function for which the locations of the two voids are in- 
dependently and uniformly distributed inside some large 
finite volume, except that the probability is set to zero 
when the distance between the void centers is less than 



2R. For this model, using the void density profile (2.4 1, 
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FIG. 15: The estimate \A1\ of the matter power spectrum 
A{k, zf for our void distribution, as a function of comoving 
wavenumber k, evaluated today at z = 0. The lower curve 
includes the correlation correction factor in square brackets in 



Eq. ( Al I, and the middle curve omits it. The upper curve is 



an approximate version of the nonlinear matter power spec- 
trum at z = obtained from the Millennium ACDM A''- 
body simulation [15], shown for comparison. The parame- 
ter values chosen were Ho = 73kms"^Mpc"\ Q.m = 0.3, 
/o = /(O) = 0.9, R = 35Mpc. 



we find for the power spectrum* 



A(fc,z) 



2a 
3^ 



1 - I2a 



,h{2kR) 



kR 



(Al) 

Here a is the void packing fraction, which is 7r/6 in our 
model, k is wavenumber, ji and j2 are spherical Bessel 
fimctions of the first kind, and /(z) is the fraction of 
the void mass in the shell, which can be computed as 
a function of redshift using the results of Sec. |II A| We 
note that this power spectrum is not an exact represen- 
tation of our void model, because in our procedure we 
first choose a direction to the source and then generate 
a density perturbation field that depends on this direc- 
tion. Thus, our procedure does not correspond exactly 
to choosing a direction randomly in a pre-existing homo- 
geneous, isotropic random process^, i.e. {Sp{x)Sp{y)) is 
not just a function of |x — y|. Homogeneity is necessary 
in order to represent the two point function in terms of 
a power spectrum. 

The power spectrum (Al) is shown in Fig. 15 both 
with and without the correction factor in square brack- 
ets that arises from the correlation between void loca- 
tions. For comparison, we also show in Fig. [15] an esti- 



* This model is not completely consistent, since the power spec- 
trum can become negative for large packing fractions. The incon- 
sistency is presumably a signal that our assumed 2- void probabil- 
ity distribution cannot be obtained starting from any symmetric 
non-overlapping n-void probability distribution. We ignore this 
inconsistency here since the correlation effects that give rise to 
the correction factor in square brackets in Eq. | |A1| | give only a 
small (< 1%) correction to (ft^) in any case. 

^ If the model were exactly homogeneous there would be a nonzero 
probability for the observer to be located inside a void. 



mate of the nonlinear power spectrum^" obtained from 
the Millennium simulation [15]. The figure shows that 
our assumed void model is in rough agreement with the 
simulation: the two power spectra agree to within a fac- 
tor ^ 2—3 at large scales, for 3 Mpc ^ ^ 30 Mpc, but 
disagree at small scales fc^^ <C IMpc, where the Millen- 
nium spectrum contains more power. This is as expected 
because our model does not attempt to model structure 
on these small scales. 

We now turn to computing the effects of lensing using 
these power spectra. From the formula (2.25) for lensing 



convergence k in weak lensing theory, it follows that for 
subhorizon modes the variance in k is [HI [47| 



dink 



9tt 



dx w(x, Xs) 



2A(fc,z)2 



k 



(A2) 

where x is comoving coordinate, Xs is the position of the 
source and w = (1 -I- z) Hqx {xg — x) /xg is the lensing 
efficiency factor. The corresponding standard deviation 
in magnitude shift Am is am — 5-\/ (k^) / In 10, from Eq. 
(3.4). We compute the integrand of the Infc integral by 



numerically integrating over redshift, for a source redshift 
of z. 



1. The result is shown in Fig. 16 



Consider first the result for our void distribution. Fig. 



16 shows that the envelope of dl^n^) /d\nk asymptotes 
to a constant at large fc, indicating a logarithmic diver- 
gence in the variance (k^). As discussed in the body of 
the paper, this divergence is an artifact of our use of dis- 
tributional density profile for each void, with a ^-function 
on the void's surface. The divergence can be regulated 
by endowing each shell with some small finite thickness 
Ar, which is approximately equivalent to truncating the 
integral over k in Eq. (A2| at fc ^ 1/Ar. Integrating Eq. 



(A2) between 10"^Mpc~^^and lO^Mpc ^ gives the result 



CTm = 0.011, which is substantially smaller than the re- 
sult CTm = 0.031 obtained from our nonlinear method in 
The agreement is improved if we inte- 
corresponding the effective cutoff 



Sec. [IV] above, 
grate up to lO^Mpc" 
lengthscale in our simulations estimated in Sec. Ill (even 
though this shell thickness lengthscale is unrealistic). In 
this case (7™ = 0.016, a factor of ~ 2 smaller than our 
simulations. The factor ~ 2 disagreement is not too sur- 
prising, since as mentioned above the derivation of Eq. 
( A2 ) requires the assumption that the density perturba- 



tion is a homogeneous isotropic random process, which 
is violated to some extent by our void model. 

It is also of interest to compute the standard deviation 
am toT the Millennium simulation spectrum. Figure 16 



We use the following fit to the Millennium power spectrum, ob- 
tained from Fig. 9 of Ref. US): A{k,zY = o(fe) (1-1- where 
the functions a and B are chosen so that A(fc)2 = 1.40889 -|- 
1.67105X - 0.11816x2 - 0.0356049x3 - 0.0367596x"' at z = 
and A(fc)2 = 0.87558 -I- 1.56132a: - 0.117482x2 - 0.0299214x3 - 
0.0383988X'* at z = 0.98, where x = logio(fc Mpc//i). This fit is 
accurate to ~ 30%. 



18 



LogioCdi^/dlnk) 




MM. Logi(j(kMpc/h) 

I'li I'i 



Da = a{S)DA, where a is the scale factor. We now define 
the magnification relative to FRW to be^^ 



n2 



(B2) 



where Da^o is the angular diameter distance computed in 
the unperturbed FRW model. Expressing the two angu- 
lar diameter distances in Eq. ( B2 ) in terms of the confor- 
mally transformed versions, the factors of a(S) cancel^"^, 
and we obtain that 



FIG. 16: The variance of the lensing convergence per unit log- 
arithmic wavenumber, d(K^) /dink, for a source at redshift 
Zs = 1, computed from the spectra shown in Fig. |15| The 
upper curve is the Millennium simulation, the lower curve is 
our void model. 



shows that the variance of the lensing convergence per 
unit logarithmic wavenumber d^K^^ /d\nk peaks at fc ~ 
100 kpc (in agreement with Sec. 10.5 of Ref. fE'). This in- 
dicates that lensing is dominated by galactic scale struc- 
tures, as claimed by Holz & Wald [2]. The total standard 
deviation^^ from all scales 10~^Mpc < < lO^Mpc is 
am = 0.044. The standard deviation from integrating 
only over the scales of voids 3Mpc < fc~^ < lO^Mpc is 
am — 0.010, a factor ~ 4 smaller; this standard deviation 



(B3) 



agrees well with our estimate ( 1.2 ) for the thick-wall void 
model. 



Appendix B: Derivation of procedure for computing 
magnification distribution 

In this appendix we describe in more detail the deriva- 
tion of our prescription for computing magnifications 
along a ray given by Eqs. (2.12) - (2.161. 

Consider an observer O and a source S. The angular 
diameter distance Da{0,S) is defined by 



D\ = SA/sn 



(Bl) 



where 5A is the proper area of the source, orthogonal to 
the direction to the observer, and 6n, is the observed solid 
angle at the observer subtended by the source. Under a 
conformal transformation of the metric, SQ is invariant 
while 6A transforms by a factor of the conformal fac- 
tor evaluated at the source. It follows that if we define 
Da to be the angular diameter distance computed in the 
conformally transformed spacetime (2.13), then we have 



where Xg is the comoving coordinate of the source. 

To compute the angular diameter distance Da{0,S) 
in the conformally transformed spacetime ( 2.13[ ), we use 
the same method that Holz & Wald used in the phys- 
ical spacetime, whose derivation we now outline in the 
context of an arbitrary spacetime. Let k — d/dx be the 
past-directed tangent vector to the null geodesic joining 
O and S, where x is affine parameter with a; = at C 
We choose vectors I, 61,62 at O so that = {k,l,eA), 
A = 1,2 form an orthonormal basis, i.e., satisfy fc^ — 
P ^ k ■ ca ~ I ■ eA = 0, k ■ I = —1, ba • eg — 6ab- 
This orthonormal basis is extended along the geodesic 
by parallel transport. 

Now let 7f{x) be an infinitesimal connecting vector that 
joins the geodesic to some nearby geodesic. The compo- 
nents of if on the orthonormal basis satisfy the geodesic 
deviation equation d^ri"/dx^ = —W^^^^kgk^-qj. More 



explicitly, expanding ff — fik + i/l 
deviation equation becomes 



rj^BA, the geodesic 



0, 



(B4a) 
(B4b) 
(B4c) 



Here dots denote derivatives with respect to x, TZ = 
Rabcdk''l''k^l'', Ua = -Rabcdk^l^k^e^, and Uab = 

Rabcdk"^ B A^'^ ^ B ' 

We are interested in a set B of rays all of which pass 
through O and which define an element of solid angle 
Sil at O. The corresponding deviation vectors 77(0) must 
vanish at O, and the initial derivatives dff/dx{0) are or- 
thogonal both to k and to the four velocity of the ob- 
server, uo- If we specialize the choice of orthonormal 
basis so that uo ■ ca 0, then it follows that = = 



This total standard deviation due to lensing computed using 
weak lensing theory and the Millennium simulation agrees well 
with that computed using other methods. For example, the cor- 
responding standard deviation for Zg = 1.5 = 0.066, which 
agrees within ~ 20% with the standard deviation of the distri- 
bution shown in Fig. 1 of Ref. (4]. 



This definition could equivalently be expressed in terms of lumi- 
nosity distances D^, since = (1 -I- z)^Da for any spacetime. 
We neglect the contribution to fi caused by the perturbation in 
the observed redshift of the source, which enters when we express 
the magnification in terms of the observed redshift. This effect 
gives a subdominant contribution to fi for subhorizon modes [481 

HS]. 
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at O, and from Eq. (B4a) we obtain that ^{x) = every- 



where. By the Unearity of the geodesic deviation equation 
it now follows that 



(B5) 



for some 2x2 matrix A\. This matrix satisfies the 



differential equation (2.141 and initial conditions (2.15) 



given in Sec. II C above, from Eq. (B4c) with v = 0. We 



define the quantity 



A{0,S) = 



dctA{xs) ' 



(B6) 



which is the so-called van Vleck determinant [5D]. One 
can show that this is invariant under rescaling of affine 
parameter, under changes of the orthonormal basis that 
preserve fc, and under interchange of O and S. 

We now define a set of angular coordinates 6^9^ that 
parameterize the solid angle measured by the observer, 
by 9^ = 9q + r]^{0)/{k ■ uo), where Oq is the direction 
to the source. The element of solid angle is then 



1 



{k ■ uoY Jb 
1 



(fc • uqY |det^(xs)| JB 



d'j^^{x,), (B7) 



where we have rewritten the integral using the Jacobian 



of the transformation ( B5 ) 



Now consider the element of area 5 A measured at the 
source S. This is defined to be the area in the rest frame 
of the source, orthogonal to the direction to the observer. 
We choose an orthonormal basis fc, l\ at 5 so that the 
four velocity is (fc + Z')/2, and decompose the connect- 
ing vector as fj = fi'k + v'l' + r( e^. Then the area 
is just &A = Jg d^T]' . Now the two orthonormal bases 

{k,l,eA) and (fc,r,e^) at S are related by some fixed 
Lorentz transformation, so we obtain 



B 



= f^' +^-i.'r>' + HABv'^D^, 



V = 



TT B lA 

Ha V 



(B8) 



for some 50(2) matrix Hab and vector D^. Since v = 
everywhere it follows that rj^ and if^ are related by 
an S0{2) transformation, which preserves ar ea, a n d so 



5 A =Ld'^ri^[xs). Combining this with Eqs. (Bl), (B6) 



and (B7) now gives for the angular diameter distance 



DA{0,Sf = 



xl(k ■ uqY 
|A(0,5)| 



(B9) 



This is independent of the normalization of the affine 
parameter and of the four- velocity of the source, but does 
depend on the four-velocity of the observer. 



We now apply the formula (B9) to a stationary ob- 
server in the perturbed Minkowski spacetime ( |2.13 ), to 
obtain the angular diameter distance Da of Eq. (B3) 
above. Specializing the affine parameter x to be the co- 
moving coordinate gives k ■ uq = 1, and then combining 



Eqs. (B3|, (B6) and (B9| gives the magnification formula 

(pjel). 



Finally, we note that in computing the matrix ^(xs), 
we follow Holz & Wald [2] in neglecting the infiuence 
of the metric perturbation on the background geodesic, 
and on the parallel transport of the orthonormal basis. 
The corresponding corrections to the angular diameter 
distance have been computed in the weak lensing limit in 
Refs. [48ll49] and are subdominant for subhorizon modes, 
that is, are suppressed by a factor of {HqR)^. 



Appendix C: Comparison with other studies of 
lensing due to voids 

Luminosity distance in the context of Swiss Cheese cos- 
mology has been studied by Clifton & Zuntz [28^ , Brouza- 
kis, Tetradis & Tzavara [23111], Szybka [3D], Valkenburg 
j31t l32i and Biswas & Notari '33] . Other studies in per- 
turbed FRW cosmologies have been done by Holz & Wald 
2J and Hui & Greene [51) . In this appendix we summa- 
rize the relevant results from this literature and compare 
with our results. 

In Clifton & Zuntz [55], the mean and standard devi- 
ation of apparent magnitude shifts are studied for red- 
shifts up to Zs ^ I in ACDM cosmology. One difference 
between their study and ours is that they model voids 
using a fully relativistic Lemaitre-Tolman-Bondi model 
with a smooth choice of density profile, whereas we use 
a simpler Newtonian model where each void consists of 
a central uniformly underdense region surrounded by a 
zero thickness shell. Fractional corrections to the New- 
tonian approximation scale as (HqR)^ ^ 0.0001 for 35 
Mpc voids, so a fully relativistic void model is not really 
necessary; our model is substantially simpler than theirs. 
A second difference between the two studies is that they 
choose a configuration of voids where the void centers 
lie along the line of sight. Due to this choice, the lens- 
ing contributions from successive voids add coherently 
instead of random walking, which significantly changes 
the magnification probability distribution. Specifically, 
for Zg — 1 and deep voids, Clifton & Zuntz obtain a 
standard deviation in modulus shift of (Jm ^ 0.01 (their 
Fig. 16), similar to our value, but they obtain a mean 
shift of {dm) ~ 0.02, a factor ~ 10 larger than ours. 
This difference arises from their lack of randomization of 
impact parameters. 

Other similar studies are those of Brouzakis, Tetradis 
& Tzavara |29J and Biswas & Notari [33] . Brouzakis et al. 
also use a fully relativistic Lemaitre-Tolman-Bondi void 
model with a smooth choice of density profile. They find 
values of standard deviation cr„ which agree to within 
~ 30% with our model; see their Fig. 5 which applies 
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to i? = 40 Mpc voids at = 1. Brouzakis et al. 
and also Biswas & Notari [33] studied the dependence of 
the magnification distribution on void sizes, source red- 
shift, and fractional underdensity in the void interior, 
and found results which agree qualitatively with ours. 
The effects of randomizing void impact parameters was 
also studied by Szybka [3D] , who found as did we that the 
dimming effect due to voids is not enough to mimic the 
effect of dark energy. The effect of shear is also studied by 
Szybka, who found its effects to be very small, in agree- 
ment with our results discussed in Sec. IIVBI above. The 
main advantage of our model compared to these studies 
is simplicity: our model allows us to explore and under- 
stand the effects of a wide range of parameter values. 

Kainulainen & Marra [THl [H] introduce a different 
technique to study lensing. While we compute the prob- 
ability distribution of magnifications by doing Monte 
Carlo simulations of ray tracing, Kainulainen & Marra 
[T5] develop a method that allows them to rapidly com- 



pute an approximate form of the entire probability distri- 
bution through a combination of numerical and analyti- 
cal techniques. However, their application of this method 
focus on the lensing due to galaxies and halos, not on the 
larger-scale structures of sheets and voids, so our study 
is not directly comparable to theirs. We note however 
that it should be possible to apply their techniques to 
compute the lensing due to voids. 

Finally, a recent paper by Lavallaz & Fairbairn |52j 
performs a similar study modeling voids as 30 Mpc 
Lemaitre-Tolman-Bondi spheres with Kostov parameter- 
ization [55] . They assume that the supernovae number 
density is proportional to the mass density inside voids 
and they study the redshift range 0.01 < z < 2.0. They 
find that if there is essentially no cut off in the lower 
range of z, the scatter in the inferred equation of state 
parameter w is about 10%, while imposing a cut off in 
the lower range of z decreases the scatter. 



[1] J. Wambsganss, R. Cen, G. Xu, and J. P. Ostriker, As- 

tron. J. Lett. 475, L81+ (1997). [21 

[2] D. E. Holz and R. M. Wald, Phys. Rev. D 58, 063501 

(1998). [22 

[3] P. Valageas, Astron. Astrophys. 354, 767 (2000). 

[4] D.E. Holz and E. V. binder, Ap.J. 631 678 (2005). [23 

[5] D. Munshi et al., Phys. Rep. 462, 67 (2008). 

[6] D. Sarkar, A. Amblard, D. E. Holz, A. Cooray, Ap. J. [24 
678, 1 (2008). [25 

[7] D. E. Holz, S. A. Hughes. ApJ 629, 15 (2005). [26 

[8] G. Ghirlanda et al. New J Phys 8, 123 (2006). [27 

[9] L. Amendola et al., Phys. Rev. Lett. 105, 121302 (2010). 
[10] C. M. Hirata, D. E. Holz and C. Cutler, Phys. Rev. D [28 
124046 (2010); C. Shang, Z. Haiman, Mon. Not. R. As- 
tron. Soc. 411, 9 (2011). [29 
[11] A. Cooray, D. E. Holz, D. Huterer, Ap. J. 637, L77 

(2006); E. Linder, JCAP 0803, 019 (2008); J. Jonsson [30 
et al., A.&A. 487, 467 (2008); A.R. Zentner, S. Bhat- [31 
tacharya, Astrophys. J. 693, 1543 (2009). [32 
[12] J. Jonsson, T. Dahlen, A. Goobar, E. MortseU, A. Reiss, [33 
J. Cosmo. Astropart. Phys. 6, 002 (2007). [34 
[13] M. Bartelmann, P. Schneider. Phys Rept 340, 291-472 

(2001). [35 
[14] Y. Wang, D. E. Holz, D. Munshi. Astrophys J 572, L15- 

L18 (2002). [36 
[15] V. Springel, S. D. M. White, A. Jenkins, C. S. Frenk, N. [37 
Yoshida, L. Gao, J. Navarro, R. Thacker, D. Croton, J. [38 
Helly, J. A. Peacock, S. Cole, P. Thomas, H. Couchman, 
A. Evrard, J. Colberg, F. Pearce. Nature 435, 629-636 [39 
(2005). 

[16] E. Lawrence et al., Ap. J. 713, 1322 (2010). 

[17] B. Li et al., arXiv: 1012. 1625 (astro-ph); S. Hilbert, J.R. [40 

Gair, L.J. King, Mon. Not. R. Astron.Soc. 412, 1023 

(2011). [41 
[18] K. Kainulainen and V. Marra, Phys. Rev. D 80, 123020 

(2009). [42 
[19] K. Kainulainen, V. Marra. Phys Rev D 83, 023009 

(2011). 

[20] R. A. Vanderveld, E. E. Flanagan, I. Wasserman. Phys 



Rev D 78, 083511 (2008). 

N. Sugiura, K. Nakao, D. Ida, N. Sakai, and H. Ishihara, 
Prog. Theor. Phys. 103, 73 (2000). 

T. Kai, H. Kozaki, K. Nakao, Y. Nambu, and C.-M. Yoo, 
Prog. Theor. Phys. 117, 229 (2007). 
N. Brouzakis, N. Tetradis, and E. Tzavara, JCAP 0702, 
013 (2007). 

B. Jain and U. Seljak, Ap.J. 484, 560 (1997). 

B. Jain et al., Ap. J. 530, 547 (2000). 
T. Hamana et al., Ap. J. 529, 56 (2000). 

T. Hamana et al., Mon. Not. R. Astron. Soc. 330, 365 
(2002). 

T. Clifton, J. Zuntz. Mon. Not. R. Astron. Soc. 400, 2185 
(2009). 

N. Brouzakis, N. Tetradis, E. Tzavara, JCAP 0804, 008 
(2008). 

S. J. Szybka. arXiv:1012.5239vl [astro-ph.CO]. 

W. Valkenburg. JCAP 0906 (2009) 010. 

W. Valkenburg. arXiv: 1104. 1082 [gr-qc]. 

T. Biswas, A. Notari, JCAP 0806, 021 (2008). 

V.Marra, E.W. Kolb, S.Matarrese, A. Riotto. Phys Rev 

D 76, 123004 (2007). 

V. Marra, E. W. Kolb, S. Matarrese. Phys Rev D 77, 
023003 (2008). 

V. Marra. arXiv:0803.3T52]>^3 [astro-ph] (2008). 

F. Hoyle, M. S. Vogeley ApJ 607 751 (2004). 

H. Mathis, S. D. M. White. Monthly Notices of the Royal 

Astronomical Society 337, Issue 4, 1193-1206 (2002). 

S. Gottlber, E. L. Lokas, A. Klypin, Y. Hoffman. Monthly 

Notices of the Royal Astronomical Society 344, Issue 

3,715-724 (2003). 

M. Plionis, S. Basilakos. Monthly Notices of the Royal 
Astronomical Society 330, Issue 2, 399-404 (2002). 
D. C. Pan, M. S. Vogeley F. Hoyle, Y.-Y Choi, C. Park. 
|arXiv:1103. 4156y2 [astro-ph.CO] 

C. Adami, A. Mazure, M. Pierre, P.G. Sprimont, C. Lib- 
brecht, F. Pacaud, N. Clerc, T. Sadibekova, J. Surdej, 
B. Altieri, P.A. Due, G. Galaz, A. Gueguen, L. Guen- 
nou, G. Hertling, O. Ilbert, J. P. LeFvre, H. Quintana, I. 



21 



Valtchanov, J.P. Willis, M. Akiyama, H. Aussel, L. Chi- 
appetti, A. Detal, B. Garilli, V. LeBrun, O. LeFvre, D. 
Maccagni, J.B. Melin, T.J. Ponman, D. Ricci, L. Tresse. 
'arXivilOlO.eigSvl [astro-ph.CO] 
[43] K. Kreckel, E. Platen, M. A. Aragn-Calvo, J. H. van 
Gorkom, R. van de Weygaert, J. M. van der Hulst, K. 
Kovac, C.-W. Yip, P. J. E. Peebles. arXiv:1008.4616vl 
[astro-ph.CO] 

[44] Margaret J. Cellar, J. P. Huchra. Science 246, 4932, 897 
- 903, (1989). 

[45] S. D. M. White, C. S. Frenk, M. Davis, G. Efstathiou. 

ApJ 313, 505 - 516 (1987). 
[46] A. Cooray, D. Huterer, D. E. Holz, Phys. Rev. Lett. 96, 



1301 (2006). 
[47] P. Valageas. A & A 356, 771 (2000). 
[48] C. Bonvin, R. Durrer, M. A. Gasparini, Phys. Rev. D 73, 

023523 (2006). 

[49] E. E. Flanagan, E. Rosenthal, I. Wasserman, Phys. Rev. 
D 79, 044032 (2009). 

[50] M. Visser, Phys. Rev. D 47, 2395 (1993). 

[51] L. Hui, P. B. Greene. Phys Rev D73, 123526 (2006). 

[52] A. de Lavallaz, M. Fairbairn, [arXiv: 1106. 161 1/ 1 [astro- 
ph.CO] 

[53] V. Rostov, |arXiv:0910.2611 



